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Recent achievements in trapping quantum gases of ultra-cold atoms have opened 
up many exciting possibilities for the experimental investigation of quantum effects in 
low-dimensional many-body systems. Among these, progress in realizing interacting 
Fermi gases promises new opportunities for studying the pronounced crossover from a 
Bardeen-Cooper-Schrieffer (BCS) superfluid to a Bose- Einstein condensate (BEC) p. 
A crossover from weakly attractive Cooper pairs to strongly repulsive molecular dimers 
can also take place in a confined quasi ID Fermi gas of atoms by using magnetic field 
induced Feshbach resonances [21 E]- In this scenario, a number of key ID integrable 
models with a variable interaction parameter have been used to describe the crossover 
0, IH El El 13 El E] • These include the integrable interacting Fermi gas [TO] . 

The ID Fermi gas can be experimentally realized [2] by tightly confining the 
atomic cloud in the radial directions and weakly confining it along the axial direction 
0, HI El El ELL]- To date, the theoretical study of the ID continuum Fermi gas has been 
mainly for the zero polarization case in the thermodynamic limit (TL). Our aim here is 
to present the results of a comprehensive study of the ID exactly solved two-component 
Fermi gas with arbitrary polarization, throughout the whole interaction range. This 
includes finite systems in the weak and strong coupling regions, and in the TL. 

The model. The Hamiltonian [10] 



describes a 5-function interacting gas of N fermions each of mass m constrained by 
periodic boundary conditions to a line of length L. The coupling constant giD can be 
written in terms of the scattering strength c = 2/<2id as gm = h 2 c/m. An effective ID 
scattering length aiD can be expressed through the 3D scattering length for fermions 
confined in a ID geometry. We use the dimensionless coupling constant 7 = mgii)/(h 2 n) 
for physical analysis. Here n = N/L is the linear density. The coupling guy can vary 
from —00 to 00 by detuning the 3D scattering length. The interaction is attractive for 
#id < and repulsive for g 1D > 0. 

The energy is given by E = ^ YljLi kj > where the quasimomenta kj satisfy the 
Bethe equations (BE) [10] 



Here j = 1, . . . , N and a = 1, . . . , M, with M the number of spin-down fermions. For 
convenience we take N even. Of particular interest are the bound states for attractive 
interaction. In this regime, the model ([Tj) has been extensively studied with N = 2M 
in the context of the BCS-BEC crossover [31 HI El El E] • For weak attraction, the system 
describes weakly bound Cooper pairs. In the strongly attractive regime the bound states 
behave like bosonic tightly bound molecular dimers. On the other hand, for repulsive 
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interaction, there are no bound states for the ground state. The spin-up and spin-down 
fermions behave like indistinguishable free fermions as c — > oo. The ground state energy 
per particle in the TL, i.e., N, M, L — * oo, N/L, M/L constant, can be derived from the 
Gaudin integral equations [10J. Here we first analyze the BE directly, i.e., for a finite 
system. 

The ground state. We observe that for weak attractive interaction a pair of fermions 
with spin-up and spin-down states forms a bound state. On the other hand, for weak 
repulsive interaction the fermion momenta separate in quasimomentum space. For 
L\c\ <C 1, it is convenient to distinguish between unpaired kj (j = 1, ... ,N — 2M) 
and paired kj (j = 1, . . . , 2M) momenta. On expanding the BE in powers including 
0(c) p2], one obtains k^ = vr(M — N — 1 + 2j)/L + 5] u) , j = 1, . . . , N/2 — M, 
kf = tt(3M - N - 1 + 2j)/L + d\ n \ j = N/2 - M + 1, . . . , N - 2M, and kf = 



7T 



1 — M + 2j+) / L + 6 fj ± v^/I, (j = 1, . . . , M), where j+ = j, if j odd and 



j + = j — 1 if j even. Therefrom it follows that the binding energy in the weakly 



attractive case is c/L for finite systems. The deviations 5 from /c| u,p ' ) are linear in c, 



with 



V = T^'t^T. (3) 



k {u) - k (p) 



r — / 



1 



^ i.(p) _ i.(p) 2 ^ b (p) - k {u) 



(4) 



where the primed sum counts each k^ only once and excludes k^. Here kj 1 ^ are 
the quasimomenta of non-interacting particles as given above. Note that the unpaired 
momenta do not interact with each other, which is consistent with the Pauli principle. 
On the other hand, for strong attractive interaction, i.e., Lc <C — 1, tightly bound states 
with quasimomentum k^ ~ Aj±|ic separate from unbound states with quasimomentum 
k (u) = n^r ^ + iuyi w]ih integerg n . = ±1> ±3 ; . . . ; - 2M - 1). In the above, 

A t = ^(l + jt + 2(jV ~ e 2A/) j for m = ±1, . . . , ±|(M - 1). In this scenario the bound 
states behave like hard-core bosons due to Fermi pressure. 

We now define the fraction a = 2M/N which characterizes the polarization of 
the system^ The ground state energy per particle in the weakly interacting regime 
follows from the above asymptotic roots of the BE as (in units of h 2 n 2 /(2m), with 
a = (1 - a) ((1 - \a) 2 + \a)) 

E a 2 . . 7T 2 n 7T 2 2 . . . 

Jj= y7 + 7(l-a)«+^ a 3 + y«o + 0(7 2 )- (5) 

The first two terms represent the binding and the interaction energy, respectively, which 
are both linear in 7. This is different from what is obtained in the TL from the Gaudin 
integral equations [H El [7], where the binding energy is oc 7 2 . The other terms in (jSJ) 
are the kinetic energies for paired and unpaired fermions. In the special case a = 0, 

I The polarization dependence of atom- hole condensation in optical lattices was discussed in Ref. [13| . 



Exact results for the ID interacting Fermi gas with arbitrary polarization 



4 



the model reduces to fully polarized free fermions without s-wave scattering. Using the 
given asymptotic expressions for the Bethe roots in the strongly attractive regime, we 
obtain (in units of h 2 n 2 /(2m)) 

- 77 2 + ^ it-a (1 - 26i) + - 2b 2 )) + 0(l/ 7 2 ), (6) 



N 4 ' 3 \16 

where ai = (1 — a) 3 , b\ = 0/(27) + 2(1 — a)/ 7 and b 2 = 20/7. Whereas the first 
term represents the binding energy, the remaining interaction energies stem from the 
interaction between pair-pair and pair-unpaired fermions. The ground state energy is 
lowest for a = 1, compared to the others for < a < 1. This is in general agreement 
with the Lieb-Mattis theorem [14| . 

Now consider the crossover to the TL. First we observe that in the TL, the 
results (J3]), (EJ) reduce to those derived from the Gaudin integral equations [U [3, EJ [7j . 
Furthermore, from ([3]), (j3J), we obtain the quasimomentum distribution functions n UiP (/c) 
per unit length in the TL, 

with A = txM/L and B = n(N + M)/L. Especially, for M = N/2, n u = 0. In Fig. [H 
n u , P (k) as calculated from the finite BE are compared to Eqs. ([7]), (jHJ). 

We now use the Gaudin integral equations to numerically compute the dependence 
of E/N on 7 for different a. The result is shown in Fig. [2j The numerical data agree with 
(jHJ) and (EJ) in the weakly interacting and strongly attracting regimes. A straightforward 
analysis of the integral equations in the strongly repulsive regime leads to 

E_ [ ^^(l-^)+0(l/7 2 ,a 2 /7), a«l 



which provides a further check on the numerics. As can be seen from ([9]), for 
strongly repulsive interaction the mixture behaves like a single component Fermi gas 
of indistinguishable particles, rather than hard-core bosons [HEllT]. 

Spin-charge separation. In the TL, the low-energy behaviour of interacting Bose and 
Fermi gases is universal in the dispersion relation and correlation behaviour [15] . Spin- 
charge separation is one of the typical Luttinger liquid behaviours, with both the spin 
and charge degrees of freedom having different propagation velocities at Fermi surfaces. 
The charge- velocity can be calculated exactly from v 2 = ^fyjr [IS] for arbitrary a, 
which yields 



v. 



^(l-46 1 ) + a 1 (l-46 2 ), 7 <-l, 



./^ + f + flo, |7l«l, 



(10) 



where the Fermi velocity Vp = him/m. 
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Figure 1. The quasimomentum distribution function (crosses) for a finite system 
(N = 50, M = 1 6) for the a) weakly repulsive and b) weakly attractive regime. The 
full lines connecting the crosses are guides to the eye. The insets are zooms of the 
distribution near the noninteracting cases (dotted lines), l/2it for unpaired fcrmions 
and l/vr for paired fcrmions. Here the full lines arc the results obtained in the TL 
(Eqs. 0, ©). 



To illustrate spin-charge separation for non-zero interaction, we show the spin and 
charge velocities in Fig. [3] as a function of 7 for a — 1. These curves are obtained 
numerically from the dressed-energy formalism [T71 [18] . For weak interaction [6] , 

V F 2 V 7T 2 



l±Za), N«l, (11) 



where the upper sign refers to v c . Furthermore, ^ = 1 — ^ = ^ for 7 — > 00 
and ^ = |(1 — I/7) for 7 — > —00 [6]. In the attractive regime, the charge excitations 
are still gapless; however, a gap opens in the spin channel [19]. This gap increases 
with increasing I7I, leading to a divergent spin velocity in the strongly attractive limit, 
with ^ = — as 7 — > —00 [6J. Some possible ways to experimentally observe the 




Figure 2. Ground state energy in units of h 2 n 2 / (2m) vs interaction strength 7 for 
different polarizations. Dashed lines compare the analytic approximations (O, 

with the numerical solution of the integral equations (solid lines) in the different 
coupling regimes. The curves shown are for a = 0, 0.2, 0.4, 0.6, 0.8, 1 (top to bottom). 



spin-charge separation in harmonically trapped atoms have been suggested [TT] . 

In conclusion, we have examined in detail the exactly solved two-component 
ID interacting Fermi gas for arbitrary polarization a in the whole interaction range. 




Figure 3. Spin and charge velocities vs interaction strength for a = 1. The straight 
dotted lines show the lowest order in 7 from (TTTj) . which is also accessible by field 
theory [TT] , 



Exact results for the ID interacting Fermi gas with arbitrary polarization 



7 



Analytic results have followed for the quasimomentum density profile, the ground state 
energy and the spin-charge velocities. These quantities reveal significant insight into 
the nature of quantum effects in the Fermi gas with arbitrary polarization. 
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